Weak and strong chaos in FPU models and beyond 
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We briefly review some of the most relevant results tliat our group obtained in the past, while 
^ I investigating the dynamics of the Fermi-Pasta-Ulam (FPU) models. A first result is the numerical 

evidence of the existence of two different kinds of transitions in the dynamics of the FPU models: i) 
a Stochasticity Threshold (ST), characterized by a value of the energy per degree of freedom below 
04 ' which the overwhelming majority of the phase space trajectories are regular (vanishing Lyapunov 

exponents). It tends to vanish as the number of degrees of freedom is increased, ii) a. Strong 
Stochasticity Threshold (SST), characterized by a value of the energy per degree of freedom at 
which a crossover appears between two different power laws of the energy dependence of the largest 
, Lyapunov exponent, which phenomenologically corresponds to the transition between weakly and 

. strongly chaotic regimes. It is stable with A''. A second result is the development of a Riemannian 

geometric theory to explain the origin of Hamiltonian chaos. The starting of this theory has been 
motivated by the inadequacy of the approach based on homoclinic intersections to explain the origin 
of chaos in systems of arbitrarily large A-, or arbitrarily far from quasi-integrability, or displaying 
a transition between weak and strong chaos. Finally, a third result stems from the search for the 
C/3 ' transition between weak and strong chaos in systems other than FPU. Actually, we found that a 

very sharp SST appears as the dynamical counterpart of a thermodynamic phase transition, which 
in turn has led, in the light of the Riemannian theory of chaos, to the development of a topological 
S ' theory of phase transitions. 

, PACS numbers: 05.45.+b; 05.20.-y 
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O In a foreword to their co-authored work reprinted in the Fermi Collected Papers [1% S. Ulam wrote: 

"...Fermi expressed often the belief that future fundamental theories in physics may involve non-linear 
operators and equations, and that it would be useful to attempt practice in the mathematics needed 
^ for the understanding of nonlinear systems. The plan was then to start with the possibly simplest such 

CN| physical model and to study the results of the calculation of its long-time behavior The motivation 

00 then was to observe the rates of mixing and thermalization with the hope that the computational 
results would provide hints for a future theory. One could venture a guess that one motive in the 
selection of problems could be traced to Fermi's early interest in the ergodic theory..." 

Actually, Fermi's early interest in ergodic theory is witnessed by his contribution to a theorem 
due to Poincare and thenceforth known as the Poincare-Fermi theorem. This asserts that neither 
analytic (Poincare) nor smooth (Fermi) integrals of motion besides the energy can survive a generic 
perturbation of an integrable system with three or more degrees of freedom, thus, in the absence of 
^ other isolating integrals of motion, any constant energy surface of these generic systems is expected 
I ' to be everywhere accessible to the phase space trajectory. At this level, no hindrance to ergodicity 
■ O seems to be possible. Whence the surprise of Fermi, Pasta and Ulam (FPU) when no apparent ten- 
^ dency to equipartition was observed in their numerical experiment whose 50th anniversary we are 
Q celebrating. Fermi himself considered what they found a "little discovery". The almost contempo- 
• • rary announcement by Kolmogorov of the starting of what would later become the celebrated KAM 
^ theorem, seemed to provide an explanation to the unexpected FPU's results. But later developments 
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of KAM theory, including optimal estimates of the iV-dependence of the perturbation thresholds and 
the Nekhoroshev theorem, revealed that this is not really an adequate framework to explain the FPU 
problem. The rich variety of the numerical phenomenology accumulated over time seemed to keep 
off "the hope that the computational results would provide hints for a future theory". In fact, "rates 
of mixing and thermalization" have a startling and complicated dependence on energy, number of 
degrees of freedom and initial conditions. Actually, any dynamical evolution of the system depends 
on the starting point in phase space and on the "landscape" of its surroundings. Thus, there can be 
a huge variety of dynamical behavior entailed by the preparation of the system in an initial condition 
out of equilibrium. As a consequence, in order to get some global information on the phase space 
structure, independently of the initial conditions, one has to look at the chaotic component of phase 
space. This way of tackling the FPU problem is very illuminating and leads to the conclusion that the 
FPU problem does not threaten the validity of statistical mechanics. Moreover, this has stimulated 
the starting of a new theory of Hamiltonian chaos. 



I. INTRODUCTION 



Few problems have been studied so extensively over the last decades as the one devised originally by E. Fermi, 
J. Pasta and S. Ulam (FPU) in 1954 Q. Their purpose was to check numerically that a generic but very simple 
non-linear many particle dynamical system would indeed behave for large times as a statistical mechanical system, 
that is it would approach equilibrium if initially not in equilibrium. In particular their purpose was to obtain the 
usual equipartition of energy over all the degrees of freedom of a system, for generic initial conditions. To their 
surprise, for the system FPU considered - a one dimensional anharmonic chain of 32 or 64 particles with fixed ends 
and in addition to harmonic, cubic (a-model) or quartic (/3-model) anharmonic forces between nearest neighbors - 
this was not observed. A variety of manifestly non-equilibrium and non-equipartition behaviors was seen, including 
quasiperiodic recurrences to the initial state. In fact, a behavior reminiscent of that of a dynamical system with 
few degrees of freedom was found, rather than the expected statistical mechanical behavior. The duration of their 
calculations varied between 10000 and 82500 computation steps. These results raised the fundamental question about 
the validity or at least the generally assumed applicability of statistical mechanics to non-linear systems of which 
the system considered by FPU seemed to be a typical example. Most of the attempts to clarify this difficulty have 
approached the problem as one in dynamical systems theory. These analyses have revealed many very interesting 
properties of the FPU system and uncovered a number of possible explanations for the resolution of the observed 
conflict with statistical mechanics. The classical explanations are: i) the survival of invariant tori in the phase space of 
a quasi-integrable system (KAM theory) 0]j the existence of Zabusky and Kruskal's solitons in a special continuum 
limit leading to the integrable Korteweg de Vries equation 0,01) Hi) the existence of an order-to-chaos transition j^. 

In this paper we will first try to exhibit the reasons why this apparently bona fide statistical mechanical system did 
not behave as such and, in particular, what in our opinion the significance of this apparent failure is for the general 
validity of statistical mechanics; then, we will review a number of results that we obtained studying the dynamics of 
the FPU models or being directly inspired by it. 

There are a number of obvious questions related to the unstatistical mechanical behavior observed by FPU, which 
all address the generic nature of the results of Fermi and collaborators: 

a) Was their time of integration long enough? 

b) Was their dynamical system of iV = 32 or 64 particles in one dimension large enough, i.e. possessing a sufficient 

number of degrees of freedom, to qualify as a statistical mechanical system? 

c) Were the recurrence phenomena (to within 3%) observed by FPU, transient or generic, i.e., possibly related to a 

Poincare recurrence time? 

The search for answers to these questions made the work of FPU very seminal, spawning many new developments 
and connections in the theory of nonlinear dynamical systems, such as the connection with continuum models based 
on the Korteweg - de Vries equation, leading to solitons |3i|, heavy breathers etc., or with few degrees of freedom 
models like the Henon-Heiles 6] and the Toda lattice 0- 

Thus, the effort to resolve the so-called FPU problem has led to enormous advances in our understanding of non 
linear dynamical systems; for a review we refer to Q. Although in our opinion the FPU problem has possibly not yet 
exhausted its power of inspiration, we believe that the FPU paradox, i.e., FPU's original question, can nowadays be 
reasonably answered along the lines we are going to describe hereafter. 
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II. STOCHASTICITY THRESHOLDS IN FPU MODELS 

For some time, it has been a commonly accepted idea that the KAM theorem provides the essential answer to 
FPU's observations, i.e., for sufficiently small nonlinearities and a class of initial conditions living on non-resonant 
tori, the FPU system behaves like an integrable system and is represented by deformed tori in phase space. With 
increasing strength of the non-linearities, a progressive chaotic behavior appears, which would ultimately lead to 
the expected approach to equilibrium and equipartition. Even though we found regular regions in phase space, the 
existing typical KAM estimates of the iV-dependence of the perturbation threshold (below which a positive measure 
of KAM tori survive), are qualitatively different from our results, indicating that the physics of the FPU model is 
quite different from what is contained in these estimates. 

Thanks to the power of modern computers, we have considerably extended the calculations performed in the past 
by various authors and we have been able to reconcile different, and sometimes contradictory, aspects of the FPU 
dynamics. 



FPU-a model 



Very interesting results have been obtained revisiting the FPU-a model Q by focusing on the development of 
chaoticity in the time evolution of the system rather than on the attainment of equipartition. For these numeri- 
cal experiments, the chosen initial conditions - single mode excitations - were the same as chosen by Fermi and 
collaborators in their original experiment. 

The model is described by the Hamiltonian 0] 
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where the particles have unit mass and unit harmonic coupling constant and the end-points are fixed (gi = qN+i = 0). 

Comparing the behaviour in time of the largest Lyapunov exponent in the FPU system with that of the same 
quantity in a suitable integrable system, it has been possible to define clearly what a trapping time in a regular region 
of phase space is and to determine numerically and unambiguously its value. The integrable system we chose for this 
comparison is the Toda lattice, from which the FPU-a model can be obtained as a third order truncation of the power 
series expansion of its potential. The Toda lattice is defined by the Hamiltonian 
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The decay pattern toward zero of X'^°'^°'{t) is undistinguishable from the decay pattern of X^^^{t) up to some time 
TT, after which X^^^{t) separates from A^°'^°(t) and converges to a finite value while \^°'^°'{t) goes to zero. This 
suggests that non-integrable motions of the FPU lattice, originated by one-mode initial excitations, enter their chaotic 
component after a transient and possibly long trapping in a regular region of phase space (some kind of "Nckhoroshcv- 
like" trapping) and that equipartition is eventually attained on a finite, albeit possibly very long, time scale. 

Fig. n shows that the trapping times tt{£,N) - so defined - for the FPU-a model at different values of both the 
energy density e and of the number of degrees of freedom N, with decreasing e first tend to increase monotonically, 
then, abruptly, display an apparently divergent behavior. 

This very steep increase of tt with decreasing e suggests the existence of, at least, a very narrow bottleneck in 
phase space, through which the system can only escape with great difhculty or, perhaps, it might not escape at all. 
The sharp variation with e of the shape of tt{s) brings about a natural definition of a threshold value of e below 
which Tt seems to diverge. 

For what concerns the behavior of the largest Lyapunov exponents, when e is smaller than the threshold value, 
A^°'*°(t) and \^^^ {t) do not show any separation, even after a very long integration time. That tt is really a trapping 
time and not, for example, a trivial effect of the numerical statistics is suggested by the fact that TT{e) ~ whereas 
A(e) ~ £3/2^ ^^^^ ^ ^-1^ Fig.ia X^^"{e,N) is reported. 

The shapes of both A (e, A^) and TT{e,N) strongly suggest the existence of a threshold value - which depends 
on A ~ of the energy density, above which the motion is chaotic and below which the trajectories appear to belong 
to a regular region of phase space. This threshold is referred to as stochasticity threshold (ST). To our knowledge, its 
direct evidence in a non-linear Hamiltonian system at A ^ 2 has been found for the first time in P| . 

Fermi and coworkers chose an initial condition well below this ST (the energy density corresponding to their initial 
condition is shown by the vertical dotted line in Fig. had they taken a ten times larger amplitude of the initial 
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FIG. 1; FPU-Q model. The trapping times TT{e,N) at different values of energy density e (i.e. at different values of the 
initial excitation amplitudes), are reported. Open squares refer to the case A'^ — 32, solid triangles refer to Af = 64, open circles 
refer to A*' = 128, respectively. The endpoints of the broken lines are lower bounds for the trapping time (the cut-off of the 
integration time is at t = 4.3 x 10*). The dotted vertical line at £ = 0.00241 corresponds to the initial excitation amplitude of 
the FPU's original paper. From Ref. 0. 
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FIG. 2: FPU-a model. The largest Lyapunov exponents X{e,N) are shown for different values of the energy density e and a 
sine wave initially excited. Symbols as in Fig. here the arrows are upper bounds for A. From Ref. 1^. 



excitation, they would have observed equipartition during the integration time they used. This appears to be the 
simple but non-trivial explanation of the lack of statistical mechanical behavior observed in the original FPU numerical 
experiment. 

In order to understand whether the ST refers to a global property of the constant energy surface or is rather a 
local property of E^;, sensitive to the initial condition, two other choices of more physically generic initial conditions, 
i.e., random positions and momenta, were considered. 

In each case a threshold energy (or equivalently energy density, since N is fixed) was found. This fact suggests 
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that the phase space undergoes some important structural change as the energy is varied: we can find regions of the 
phase space where ordered trajectories are observed, regions where there is coexistence of order and chaos and regions 
where chaos is fully developed. 

An important question is whether the stochasticity threshold is stable or unstable with respect to N . Unambiguous 
information about this point is provided by the Lyapunov exponents A(e, N) computed at different N, always starting 
with random initial conditions (Fig.EJ. 
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FIG. 3: FPU-Qf model. The largest Lyapunov exponents Xie, N) are plotted vs. the energy density e, for different values of A''. 
Random initial conditions are chosen. Star-like polygons refer to N = 8, crosses to A'^ = 16, asterisks to A'' = 32 and star-like 
squares to A'' = 64, respectively. The arrows have the same meaning as in Fig. |5| From Ref. |^. 

At large e, there is a tendency of all the sets of points to join, while they tend to separate at small e: the larger 
N, the smaller the energy density at which the separation occurs. The "critical" energy density Sc at which the 
separation occurs shows the A^-dependence £c{N) cx l/N"^ . 

A qualitative agreement about the vanishing with N of the critical energy to get chaos is reported in a recent paper 
on the FPU-a model ^3 • The question of how to explain the existence and the 1 /A^^ dependence of the stochasticity 
threshold remains open. 

We thus see that revisiting the FPU-a model led to the observation of some very interesting phenomena: the 
apparent existence of regular regions in the phase space of a non-integrable Hamiltonian system with many degrees 
of freedom at large values of the anharmonic energy (even very large if compared with what they should be according 
to the KAM theory), and the existence of almost regular regions of phase space where the trajectories are trapped 
during long but finite times. The behavior of the largest Lyapunov exponent suggests that the sudden escape from 
the regular region might occur as if the trajectory would eventually find a 'hole' in its boundary. 

Moreover, the coexistence of regular regions of the phase space and of a large chaotic "sea" reconciles different and 
sometimes apparently contradictory aspects of the FPU dynamics found in the past. The lack of equipartition in the 
original FPU experiment is not representative of a global property of phase space: apparently regular, soliton-like 
structures, similar to those of Zabusky and Kruskal, have a very long, possibly infinite, life-time below the stochasticity 
threshold, whereas, above the same threshold, they have only a finite life-time 9]. 

The threshold energy density for the onset of chaos shows a clear tendency to vanish at an increasing number of 
degrees of freedom 1/N'^), so that strong evidence has been found in support of the point of view that the so-called 
"FPU-problem" does not invalidate the (generic) approach to equilibrium and the validity of equilibrium statistical 
mechanics. On the other hand the existence of long living initial states far from equilibrium, may well have interesting, 
non trivial physical applications. 
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B. FPU-/? model 



The FPU-/3 model is described by the Hamiltonian Q 



N r 



k=l 



(3) 



where the particles have unit mass and unit harmonic coupling constant and the end-points are fixed [qi = qn+i = 0); 
for this model also periodic boundary conditions have been considered (qi = qpf^i). 

The approach to equilibrium of the FPU-/3 model was studied extensively for various classes of initial conditions by 
Kantz et al. |0| and recently by De Luca et al. who extended and improved earlier computations of ours [T3.ll4j . 

A very detailed picture has emerged from these works, as to the behavior of the FPU-/? model in its dependence 
on non-equilibrium initial conditions as well as in the role played by low frequency and high frequency mode-mode 
couplings [l5( during its time evolution. 

Several years ago, we introduced a time dependent spectral entropy S{t) — — Wi{t) log Wi{t), where Wi{t) — 
Ei{t) / '^1^ Ek(t) is the normalized energy content of the i-th harmonic normal mode, defined so as to detect energy 
equipartition (when it attains its maximum value) and to measure the time needed to reach it. By means of this 
spectral entropy, we investigated in Refs. [TgL Ii7| the relationship between equipartition times, measured through 
the time relaxation patterns of this spectral entropy, and the chaotic properties of the dynamics in nonlinear large 
Hamiltonian systems. For the FPU-/3 model, we have put in evidence that, at different initial conditions and at long 
times, the spectral entropy always relaxes toward its maximum value signaling equipartition, however, depending on 
the value of the total energy density, the relaxation occurs with quite different modalities. The relaxation time is 
approximately constant for energy densities greater than some threshold value Ec, but it steeply grows by decreasing the 
energy density below this threshold. Moreover, the largest Lyapunov exponent shows a crossover in its e-dependence 
corresponding to this threshold value. We interpret this phenomenological result as the (smooth) transition - at Sc 
- between two different regimes of chaoticity, weak chaos and strong chaos, thus we called this transition the Strong 
Stochasticity Threshold (SST) O. Weak and strong chaos are qualitative terms to denote slow and fast phase 
space mixing respectively. In Refs. |l6l Il7j we resorted to a random matrix model for the tangent dynamics to try 
to make more precise and quantitative the definitions of weak and strong chaos. At least in a limited high energy 
range of values, the random matrix model predicts the numerically observed scaling A(e) ^ e^/^ (this law changes 
to A(£) - £1/4 at very high energy density, however this is not explained by the random matrix model, the reason 
is that a free parameter, a time-scale of unknown e dependence, enters the random matrix model. This time-scale 
is arbitrarily assumed constant). Thus we say that chaos is strong in the energy density range where A(e) ~ e^/'^, 
because the random matrix model assumes that the dynamics looks as a random uncorrelated process (if sampled with 
the just mentioned unknown time scale). At low energy density, the e-scaling of A is found to be steeper, A(e) ^ e^, 
so that A fastly decreases as e is lowered and is much smaller than it should be if the high energy random matrix 
prediction could be extrapolated down to low energy values. For this reason we say that here chaos is weak. Figure 
ill shows A(e,iV). 

The SST is independent of the initial conditions so it has to be ascribed to some change in the global properties 
of the phase space, for this reason it has to have major consequences on the dynamics. An interesting explanation 
based on a model for phase space diffusion is given in Rcf. Ts^ . 

The SST has been found to be correlated with changes in the transient non- equilibrium behavior (e.g., relaxation 
to equipartition) IT^ llTl IT^ . and has been found to be also correlated with stationary non- equilibrium phenomena 
like heat conduction[20j. The SST is found to be independent of the number of degrees of freedom, which makes it of 
prospective relevance for equilibrium statistical mechanics. Among the model dependent consequences of the existence 
of the SST, it is worth mentioning that in the FPU-/? model, at £ < ef^'^ high-frequency excitations yield longer 
relaxation times with respect to low frequencies. This is in agreement with the common belief that high-frequencies 
have the tendency to freeze; at £ > £f the situation is reversed. High frequency excitations yield a quicker relaxation 
than low frequencies p!d| . 

It is remarkable that the existence of the SST is not only a characteristic of the FPU-/? model. In fact, it has 
been detected in the follo wing one dimensional lattices: with diatomic Toda interactions (i.e., with alternating masses 
that break integrability) l2l||: with single- well cf)^ interactions with smoothed Coulomb interactions '2lj; with 
Lennard- Jones interactions |2lj ; in an isotropic Heisenberg spin chain [22| ; in a coupled rotators chain which displays 
two thresholds separating two regions of weak chaos (occurring at low and high energies) from an intermediate region 
of strong chaos [23l |23 | ; in a "mean-field" XY chain ^5^ and in homopolymeric chains ^2d\ . It has been also detected 
in two and three dimensional lattices, with two-wells cp^ interactions 27, 28J, with XY Heisenberg interactions ,29,.3Cj |. 
Therefore the SST seems to be a generic property of Hamiltonian systems with many degrees of freedom. 
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FIG. 4: FPU-/3 model. Largest Lyapunov exponents Ai vs. energy density e at = 128 and at different initial conditions: 
random at equipartition (circles), wave packets at different average wave numbers (squares, triangles and asterisks). From 
Ref. 0. 



C. FPU-(a + /3) model 



In the FPU-a model, the existence of a stochasticity threshold (ST) at an energy density below which the dynamics 
is regular has been observed. In the FPU-/? model, a strong stochasticity threshold (SST) above which the dynamics 
is strongly chaotic has been found. By combining these two models into the FPU- (a -I- /3) model, it is possible to 
observe the coexistence of both the ST and the SST. This model has been studied recently in [sj- It is described by 
the Hamiltonian 



N 



k=l 



\pk + \{ik+i - qkf + ^(ft+i - qkf + ^(gfc+i - qkf 



(4) 



where the particles have unit mass and a unit harmonic coupling constant and the end-points are fixed [qi — qN+i = 0). 
This model Hamiltonian, with the choice of a = 0.25 and j3 — |a^, is a fourth-order expansion of the Toda model 
(|^ . Consequently, its potential function is very close to interatomic potentials of the Morse or Lennard- Jones type 
in solids, provided that a suitably restricted energy density range is considered. Random initial conditions have 
been chosen. The results of the computation of the largest Lyapunov exponents at different energy densities and for 
different values of N are shown in Fig. El The patterns of A(e, N), therein reported, display some remarkable features. 
For small values of the energy density, there is a sudden drop of A which, in close analogy with Rcf . 9] , allows us 
to define an ST below which we can assume that the overwhelming majority of the trajectories in phase space are 
regular. This ST moves to smaller and smaller values of e as is increased. 

Around e ~ 0.8, a "knee" is observed in the pattern A(e, A^) (Fig. [SJ, due to a crossover between two power law 
behaviors, ~ at small e and ^ eX large e, where the latter has been attributed to the existence of an SST 
[l^[l7l| . This crossover is the signature of the transition from weak to strong chaos, as already discussed in [l6lfl7| . 



III. RIEMANNIAN GEOMETRY OF CHAOS IN THE FPU-/3 MODEL 



In this Section we sketch how we have analytically computed, in the limit of arbitrarily large N , the largest 
Lyapunov exponent A as a function of the energy density e for the FPU-/? model. The excellent agreement of the 
analytic outcome with the numerical results for A(e) provides a preliminary understanding of the transition between 
weak and strong chaos (SST), and strongly supports the general validity of the proposed explanation of the origin of 
Hamiltonian chaos. 
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FIG. 5; FPU-(a + /3) model. The largest Lyapunov exponents \{e, N) are shown for different values of the energy density e for 
various values of A''. Starlike squares refer to A'' = 8, asterisks to A*' = 16, open squares to A*' = 32, open triangles to A' = 64, 
open circles to A' = 128, starlike polygons to A' = 512 and crosses to A^ = 1024, respectively. Full squares refer to A^ = 32 
and excitation amplitudes A ranging from 5 to 50. Solid lines are the asymptotic scalings e'^ and at low and high energy 
density, respectively. From Ref. [31| . 



For generic non-integrable Hamiltonian systems, when the number of degrees of freedom is large, which in practice 
means already a few hundred, the whole phase space is filled by chaotic trajectories, at least at physically meaningful 
values of the energy density. Therefore, any framework of analytic description of the dynamics has to cope with chaos. 
However, even the basic question about the origin of chaos itself, in many degrees of freedom Hamiltonian systems 
seems to lack an answer. For example, all the theoretical machinery of Classical Perturbation Theory (CPT) is of 
little use if we want to deal with chaos, and so does the traditional explanation of its origin based on homoclinic 
intersections p^ . 

Until a few years ago, the "only game in town" , which seemed of potential interest to treat chaos at large N, was an 
attempt by Krylov 32] at explaining the origin of phase space mixing as a consequence of negative scalar curvature of 
suitable Riemannian manifolds whose geodesies coincide with the solutions of Newton equations of motion. Krylov's 
idea was to take advantage of some mathematical results about the stability properties of geodesies on negatively 
curved Riemannian manifolds. These results are associated with the names of Hadamard [s^l, Hedlund 34] and Hopf 
[ssl ]. Since Krylov's, other attempts have been done along the same line of thought (see e.g. the discussion in 36]), 
but none of them appeared very useful. 

More recently, we have reconsidered the Riemannian geometric approach and, with the aid of numerical simulations 
on the FPU-/3 model, we have discovered why the previous attempts failed: the dominant mechanism for chaotic 
instability in physically relevant geodesies flows is parametric instability due to curvature variations along the geodesies, 
and not necessarily geodesic flows on negatively curved manifolds |24l l36l l37l l38l l39l l40|l. On this basis, we have 
started the formulation of a Riemannian theory of Hamiltonian chaos which applies to dynamical systems described 
by a standard Lagrangian function 

L{q,q) ^^a.kq'q'' -V{q) , (5) 

where aik is the kinetic energy matrix (aik = 6ik for the usual form of the kinetic energy) or, equivalently, by the 
Hamiltonian H{p, q) — ^a^'^piPk + V{q), where the momenta are given by pi — ai^q^ . From Maupertuis' least action 
principle for asynchronous isoenergetic varied paths ^{t) with fixed endpoints 

5 ( 2Wiq, q)dt^5 I {2[E - V{q)]a,kq'ci'Y'^dt = , (6) 

where W is the kinetic energy, the equations of motion follow. Equation JHJ is equivalent to the extremization of the 
length-integral J ds where ds is ds^ = gik{q)dq^dq'^ — 2[E — V{q)]aikdq'^dq'' . In other words, mechanical trajectories 
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are geodesies of the configuration space endowed with a proper Riemannian manifold structure described by the metric 
tensor 

g^k{q) = 2[E-V{q)]a^k . (7) 

This is known as Jacobi metric and is defined in the region of the configuration space where V{q). In local coordinates, 
the geodesic equations on a Riemannian manifold are given by 



^9* ; dq^ dq 



k 



(8) 



ds d 



s 



where s is the proper time and r*^, are the Christoffel coefficients of the Levi-Civita connection associated with gik- 
By direct computation, using = {E - V{q))5ik, T"-.^ = ^5"'^{djW5k,n + dkWS^j - d^WSjk) and ds'^ = 2W^dt^, 
it can be easily verified that the geodesic equations yield 

dt^ dq^ ' ^ ' 

i.e. Newton's equations associated to the Lagrangian (jsj. 

Among other Riemannian geometrizations of Newtonian dynamics, a very interesting one is defined in an enlarged 
configuration spacetime M x B? , with local coordinates g^, . . . , g', . . . , q'^ ,q^^^), endowed with a non-degenerate 
pseudo-Riemannian metric whose arc-length is 

ds^ = g^,l, dq'^dq" = dq'dq^ - 2V{q){dq^f + 2 dq^dq^+^ , (10) 

called Eisenhart metric. The natural motions are obtained as the canonical projection on the configuration space-time 
of those geodesies for which the arclength is positive-definite and given by ds^ = (const)^ dt^. A way of measuring 
of how much a Riemannian manifold deviates from being a Euclidean manifold is provided by the degree of non- 
commutativity of the covariant derivatives which is properly defined by the Riemann-Christoffel curvature tensor 
R{X,Y) = VxVy — VyVx, where V is the Levi-Civita connection, and X,Y are tangent vectors There are 
two relevant curvature scalars: the Ricci curvature Kr in a given direction v, and the scalar curvature TZ (see [s^'l. 

There is an important relation between the curvature of a manifold and the stability of its geodesies. In fact, 
the evolution of a vector field J, called geodesic separation vector, is completely determined by the curvature tensor 
according to the equation 

V2j(s) + i?[J(s),7(s)]7(s) = 0. (11) 

This is the Jacobi - Levi-Civita equation, where V-y is the covariant derivative in the direction of the velocity vector 
7 = w. J contains the whole information on the stability - or instability - of any given reference geodesic 7(5) because 
it locally measures the distance from 7(3) of any given geodesic close to 7(5). 

Since the Jacobi equation Hll|) relates the stability of the geodesies of a manifold to its curvature, the Jacobi 
equation links stability and instability (chaos) of the dynamics with the curvature of the "mechanical" manifold [43j , 
if the metric is associated with a physical system. 

In the particular case of isotropic - or constant curvature - manifolds, Eq. Hll() becomes very simple: choosing a 
geodesic frame, i.e., a reference frame transported parallel along a reference geodesic, the Jacobi equation is written 
as 



— +KJ = 0, (12) 



ds 

where K = Kji/{N — 1) = TZ/N{N — 1), which has either bounded oscillating solutions || J|| cx cos{^fK s) or ex- 
ponentially unstable solutions |j J|| cx expW^K s) according to the sign of the curvature and thus of the constant 
K. 

As long as the curvatures are negative, the geodesic flow is unstable even if the manifold is no longer isotropic, 
and the instability exponent is greater than or equal to (— vaayLM^K))^^'^ ■ Geodesic flows on compact manifolds with 
everywhere negative curvature were studied for the first time in the classic works by Hadamard, Hedlund and Hopf 
[33> .34 .35j and many results were then established by Anosov 01 > among them the fact that such systems are ergodic 
and mixing. 

Equation H12(l is valid only when K is constant. Nevertheless in the case in which dimM — 2 (surfaces), the Jacobi 
equation - again written in a geodesic reference frame for the sake of simplicity - takes a form very close to that of 
isotropic manifolds, 



1 

^ 2 



2 , ^7^(s)J = 0, (13) 
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where TZ{s) denotes the scalar curvature of the manifold along the geodesic 7(3). The solutions of Eq. IjlSfl may 
exhibit an exponentially growing envelope even if the curvature TZ{s) is everywhere positive but non constant. For 
example, in the case of the celebrated Henon-Heiles model 0, the scalar curvature TZ, computed with the Jacobi 
metric, is always positive despite the existence of fully developed chaos above some threshold energy 's^. As a matter 
of fact, the generic condition of physically relevant systems (like coupled anharmonic oscillators on d-dimensional 
lattices) is that Ricci and scalar curvatures of the mechanical manifolds are neither constant nor everywhere negative, 
and the straightforward approach based on Eq. H12|l does not apply. 

The key point is to realize that negative curvatures are not necessary to generate chaos, while the generic non 
constancy of the curvature of mechanical manifolds (in the absence of very "exotic" hidden symmetries jS|), triggers 
parametric instability of the geodesies. Thus the exponential growth of the solutions of the stability equation H13|l . 
that is chaos, even if no negative curvature is "felt" by the geodesies. 



A. A geometric formula for the Lyapunov exponent 

In the large N case, with some simplifying assumptions, mainly that the mechanical manifolds are quasi-isotropic 
[23, 113, it is possible to derive an effective scalar stability equation resembhng Eq. l(T^ . where the role of TZ{s) is 
played by a random process, so that an analytic estimate of the largest Lyapunov exponent can be worked out. This 
effective equation is independent of the knowledge of the dynamics and has the form 



+ n{s)^^ 



(14) 



where denotes any of the components of J because now all of them obey the same effective equation of motion, and 
the squared frequency Q{s) is a gaussian random process 



n{s) = (kn) 



,1/2 



vis) 



(15) 



where ka — Kr/{N ~ 1) and {S'^kR}^ is a shorthand for j^zti {^^ K fj,; the averages (•)^ are microcanonical averages; 
77(5) is a gaussian random process with zero mean and unit variance. Our estimate for the (largest) Lyapunov exponent 
A is then given by the growth-rate of |l(?/',-0)(i)|P according to the definition 



A = lim — log : 

t^oo2t V'^(0) + V2(0) 



(16) 



The ratio (V'^ft) + V'^(O)/('0^(O) + ri^)) is computed by means of a technique developed by Van Kampen, sum- 
marized in Ref. |24|. where the following expression for A has been derived 



A 




(17) 



where fip = (fcfl)/^; fn = {S'^kfi}fj = ■^[{Kj^)fj, — (i^i?)^] and r is a time scale expressed in terms of fJo and an. The 
quantities Oq, an and r can be computed as static, i.e. microcanonical averages. Therefore Eq. (|17|l gives an analytic, 
though approximate, formula for the largest Lyapunov exponent independently of the numerical integration of the 
dynamics and of the tangent dynamics. 

A completely analytical computation of A(£) has been performed - in the thermodynamic limit - for the FPU-/3 
model (such a result first appeared in then it was refined in and for other models. We report in Fig.|Bl 

the result for the FPU case: the agreement is strikingly good. The analytic values of A agree with the numerical ones 
with errors of a few percent in a range of six orders of magnitude both in e = E/N and A, and no use of adjustable 
parameters has been made. A preliminary explanation of the existence of the SST proceeds as follows. At low e, 
the amplitude of the curvature fluctuations an is much smaller than the average curvature f2oj thus the mechanical 
manifolds are not very different from constant curvature manifolds, so that the geodesic flow has many of the features 
that it would have if it lived on a strictly constant curvature (equal to the average curvature) manifold, and, loosely 
speaking, a slow phase space filling through tortuous paths will take place: chaos is weak. Conversely, when an flo, 
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FIG. 6: FPU-/3 model. Lyapunov exponent A vs. energy density e with /3 = 0.1. The continuous line is the theoretical 
computation according to Eq. 11711 ._ while the circles and squares are the results of numerical simulations with A'' respectively 
equal to 256 and 2000. From Ref. 



we can imagine that no similarity at all will exist between the chaotic geodesic flow and its integrable counterpart 
living on a constant curvature (equal to the average curvature) manifold. As a consequence the geodesic flow can 
quickly diffuse in any direction in phase space: chaos is strong. 

Other systems for which good results have been obtained are: a one-dimensional chain of coupled rotators p^ . 
two and three dimensional classical XY Heisenberg models ^3 j two and three dimensional classical lattice (p'^ models 
|27Ll28l |. "mean-field" XY model [S^, though some adjustments are necessary in these cases. 

An important remark is in order. The geometrical theory of chaos aims at explaining what is the origin of chaos 
in Hamiltonian systems, and not at providing a recipe for the computation of Lyapunov exponents. The impressive 
success of the theory in analytically computing Lyapunov exponents for the FPU model, means that we have actually 
found the right conceptual framework and the right explanation for the existence of Hamiltonian chaos and warrants 
that any effort to further develop the theory is worthwhile. 

One has to keep in mind that the above given analytic formula for the largest Lyapunov exponent has a limited 
validity domain: that of the fundamental assumption of quasi-isotropy of the mechanical manifolds. The next step 
will be that of relaxing the assumption of quasi-isotropy by letting in nontrivial topology of configuration space. 



IV. HAMILTONIAN DYNAMICS, PHASE TRANSITIONS AND TOPOLOGY 

Though the content of this Section could appear to be somewhat far from the initial FPU problem, we have 
nonetheless sketched it in order to remark how fertile, inspiring and far reaching a systematic investigation of the 
(once) surprising behavior of the FPU dynamics has been. 

The macroscopic properties of large- Hamiltonian systems can be understood by means of the traditional methods 
of statistical mechanics. The origin of the possibility of describing Hamiltonian systems via equilibrium statistical 
mechanics are the chaotic properties underlying the dynamics. 

Above, we have observed that the crossover in the e-dependence of A phenomenologically corresponds to a transition 
between weak and strong chaos (SST), or slow and fast mixing respectively. Thus we have surmised that this transition 
has to be the consequence of some "structural" change occurring in configuration space, and thus also in phase space. 
This dynamical (mild) transition has been observed, we said above, in many other systems besides FPU. Then, some 
natural questions arise: could some kind of dynamical transition between weak and strong chaos (possibly sharper 
than the SST found in FPU models) be the microscopic counterpart of a thermodynamic phase transition? and if 
this was the case, what kind of difference in the A(e) pattern would discriminate between the presence or absence of 
a phase transition? and could we make a more precise statement about the kind of "structural" change to occur in 
configuration space when the SST corresponds to a phase transition and when it does not? 
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During the last years, after an earlier attempt in Ref. |29| where the classical XY model in two dimensions was 
considered, and its largest Lyapunov exponent was found to display some indication of the transition temperature of 
the Kosterlitz-Thouless phase transition, there has been a renewed interest in the study of the behaviour of L yap unov 
exponents in systems undergoing phase transitions, and a number of papers have appeared: see 123. l3Clj| and 

liSi, iM, iE, M, H, 52 , 53, 15 J . 

Two systems have received considerable attention in this framework; the so-called lattice ip'^ model, and the mean- 
field XY model. The lattice (p'^ model is described by the Hamiltonian 



H = Y.ip^^ + - ^^^^ + . 



2 



(18) 



where the pi are momenta conjugated to the ipi, real valued scalar variables defined on the sites of a d-dimensional 
lattice; and u are positive parameters, and the brackets (i, j) stand for nearest-neighbors. This model has a phase 
transition at a finite temperature provided that > 1. 

The mean- field XY model [s^ describes a system of N equally coupled planar classical rotators. It is defined by 
the Hamiltonian 



1 J ^ ^ 

= ^ 2^' + ^ J2 [1 - -fj)]-hY, cos ^^ . (19) 

i — l i=l 

Here pi € [0, 2tt] is the rotation angle of the i-th rotator. Defining at each site i a classical spin vector = 
(cos (^i, sin </5i) the model describes a planar (XY) Heisenberg system with interactions of equal strength among all 
the spins. The equilibrium statistical mechanics of this system is exactly described, in the thermodynamic limit, by 
mean- field theory '56']. In the limit /i — > 0, this system has a continuous phase transition. 

Through standard methods of molecular dynamics, thermodynamical observables have been computed and found 
to be in agreement with statistical mechanical predictions. The energy density (e = E/N) dependence of the largest 
Lyapunov exponent numerically found in the (p'^ model - reported in Fig. 0- shows a pattern similar to that found 
in the FPU model but now the mild transition between weak and strong chaos is replaced by an abrupt transition, a 
sharp SST: a "cuspy" point in A(£) shows up which corresponds to the critical energy locating the phase transition. 
Also the Lyapunov exponent of the mean field XY model, obtained through an analytic estimate worked out in the 
limit N ^ oo 25] by means of the above discussed geometrical theory of chaos, sharply signals the phase transition 
(see Fig.lHl). 

In both cases, it is evident that the £-pattern of the largest Lyapunov exponent clearly signals the presence of a 
phase transition; the same happens for all the other models studied in the above mentioned references. 

Then, coming to the other questions, as Lyapunov exponents are tightly related with the geometry of the mechanical 
manifolds in configuration space (as well as in phase space), we have been led to conjecture that in presence of a phase 
transition we have to go to the deeper level of topology of these manifolds to find an adequate explanation [s^ . If this 
is actually the case, we are confronted with a possible - at least conceptual ~ deepening of our understanding of the 
origin of phase transitions. In fact, the topological properties of configuration space submanifolds, mainly equipotential 
hypersurfaces — V''^{v) = {q G R^\V{q) — v} or the regions bounded by them My = {q G R'^\V{q) < v}, are 
already determined when the microscopic potential V is assigned and are completely independent of the statistical 
measures. The appearance of singularities in the thermodynamic observables could then be the effect of a suitable 
topological transition in configuration space. Several results strongly support this Topological Hypothesis and suggest 
that a phase transition might well be the consequence of an abrupt transition between different rates of change in the 
topology above and below the critical point. More details can be found in the review paper [s^ and in the subsequent 
papers: 0| where the topology of the is analytically studied for the mean-field XY model; [5^ [g^I where the 
top olog y of the is analytically studied for a trigonometric model undergoing also a first-order phase transition; 
[SOIIM where an analytic relationship between topology and thermodynamic entropy is given among other results; 
|62l l63j | where a preliminary account of a general theorem on topology and phase transitions is given. 



V. CONCLUDING REMARKS 



By chance, Fermi, Pasta and Ulam chose the initial condition, in their numerical experiment on the a-model, below 
the threshold energy of a transition between regular and chaotic motions. With a stronger initial excitation, no 
"FPU problem" would have arisen because equipartition of energy could have been observed even with the rather 
short integration time that the authors could afford 50 years ago. Apparently the observed phenomenology cannot 
be explained by the existing formulations of the KAM theory, both because of the large degree of anharmonicity 
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FIG. 7: Lyapunov exponent A vs. the energy per particle e, numerically computed for the two-dimensional 0(1) tp"^ model, 
with A'^ = 100 (solid circles), A'^ = 400 (open circles), A'^ = 900 (solid triangles), and A'^ = 2500 (open triangles). The critical 
energy is marked by a vertical dotted line; the dashed line is the power law e^. From Ref. [2^ . 
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FIG. 8: Mean field XY model: analytic expression for the Lyapunov exponent (solid curve). The curves above the transition 
are finite- A'^ results for A' = 80 (upper dashed line) and TV = 200 (lower dashed line): here A oc N~^^^ . From Ref. |25|l . 



(nonintegrability) at which the ST occurs, and because of its slow vanishing at increasing N . Since this threshold 
energy goes to zero as the number of degrees of freedom is increased, the FPU problem is not a true problem 
for equilibrium statistical mechanics. Nevertheless, the existence of possibly long-living transient noncquilibrium 
phenomena draws attention to the relevance of dynamics, initial conditions and observational time scales in order to 
assess whether dynamics can be replaced by statistics or not. Because of the cubic potential, the a-model is unstable 
above an upper bound in energy density, so FPU considered also the so-called /3-model which is well defined at any 
energy. However, in the /3-model it is hard to detect the ST because it seems to occur at a very low energy density, 
where the convergence of the largest Lypaunov exponent requires huge computational times. On the other hand, the 
/3-model displays another and much more interesting chaotic transition, that we called SST, which is a transition 
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between weak and strong chaos. Strong chaos is related with fast phase mixing and fast thcrmahzation of an out of 
equihbrium initial condition. Weak chaos is associated with a sudden increase of relaxation times of noncquilibrium 
initial conditions when the energy density is smaller than a threshold value (which corresponds to the SST). At 
sufficiently low energy density the thermalization can be so slow that the system can give the wrong impression to 
recur ad infinitum, if the observational time is not long enough. The study of the a + /3-model, which provides a 
good approximation of interatomic interaction potentials of the Morse or Lennard-Jones type, displays both the ST 
and the SST. However, only the SST is stable with A'^ and can thus be relevant for equilibrium and nonequilibrium 
statistical properties of a large class of classical many-body systems. In fact, this kind of transition seems a common 
property of many degrees of freedom Hamiltonian systems. 

The systematic investigation of the chaotic properties of FPU models - being a heavy numerical task - has become 
possible only rather recently, with the advent of modern powerful computers. The results so obtained demanded a 
satisfactory and constructive explanation of the origin of Hamiltonian chaos as well as for the reason of the transition 
between weak and strong chaos. Motivated by the need of understanding chaos in FPU models, we have started a new 
and successful theory of Hamiltonian chaos which resorts to basic concepts and methods of Riemannian geometry. 
Later on, all these findings have suggested to look at phase transition phenomena from a new point of view which, 
eventually, has inspired the development of a new theoretical approach to them, based on topological concepts. 
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